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Abstract 

We discuss the solution of the Mott transition problem in a fully frus- 
trated lattice with a semicircular density of states in the limit of infinite 
dimensions from the point of view of a Landau free energy functional. 
This approach provides a simple relation between the free energy of the 
lattice model and that of its local description in terms of an impurity 
model. The character of the Mott transition in infinite dimensions, (as 
reviewed by Georges Kotliar Krauth and Rozenberg, RMP 68, 1996, 13) 
follows simply from the form of the free energy functional and the physics 
of quantum impurity models. At zero temperature, below a critical value 
of the interaction U, a Mott insulator with a finite gap in the one parti- 
cle spectrum, becomes unstable to the formation of a narrow band near 
the Fermi energy. Using the insights provided by the Landau approach 
we answer questions raised about the dynamical mean field solution of 
the Mott transition problem, and comment on its applicability to three 
dimensional transition metal oxides. 

1 Introduction 

The idea of understanding lattice models of correlated electrons from a local 
perspective is a very intuitive one, and has been used repeatedly in many body 
physics over the years. 

A well known example is the heavy fermion problem, where a broad band 
of conduction electrons interacts with more localized f electrons, via a magnetic 
Kondo exchange interaction. In the early days of the heavy fermion problem, 
a great deal of understanding was obtained by considering the screening of an 
isolated spin by a sea of conduction electrons, i.e. studying the single site Kondo 
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effect, and then regarding the Kondo lattice as a collection of Kondo impurity 
models. 

This view had some successes in explaining the origin of a non perturbative 
energy scale, the Kondo temperature, where the properties of the system change 
dramatically (for instance the susceptibility crosses over from Curie to Pauli 
behavior ). 

In an early paper, however, Nozieres jlj] pointed out, that the physics of 
the lattice problem is more complex than the single impurity problem since at 
least in the limit of low density of conduction electrons, there are not enough 
itinerant electron spins, to screen all the impurity spins in the lattice. In this 
case one cannot regard the lattice as a collection of single Kondo impurities. 
This issue is now known as Nozieres' exhaustion problem. This is perhaps one 
of the earliest warnings that single impurity thinking can be misleading if it is 
applied uncritically to lattice problems involving a correlated degree of freedom 
at each site. 

In the context of transition metal oxide physics, an impurity view of the 
d-electron spectral function was put forward by Zaanen, Sawatzky and Allen Q 
and by Fujimori, Minami and Sugano|j|], and led to a qualitative description of 
the spectra in these systems. In the light of the modern developments of the 
dynamical mean field theory, we would regard the early applications of impurity 
views to the physics of f and d electron systems as local (but not self consistent) 
impurity approximations to lattice models. 

The last ten years have witnessed dramatic progress in the theory of corre- 
lated electron systems. The modern developments of a dynamical mean field 
theory Q and its implementation via mappings onto impurity models [El , now 
allows us to use impurity models supplemented by a self consistency condition 
to study lattice models. The results are exact in the well defined limit of infinite 
lattice coordination [0. 

We are now in a much better position to gauge the reliability of the argu- 
ments based on the Local Impurity Approximation by studying lattice models 
in the limit of large lattice coordination using Local Impurity Self- Consistent 
Approximations. If the self consistency condition does not play an important 
role, naive impurity based arguments are reliable. Since the dynamical mean 
field theory is exact in the limit of large lattice coordination, we can also under- 
stand which physical elements are absent in this limit (a most notable example 
is the feedback of the magnetic correlations on the single particle properties) 
and assess in which physical circumstances it provides reliable guidance to the 
physics of three dimensional real materials. 

This paper is devoted to the problem of the Mott transition, i.e. the interac- 
tion driven metal insulator transition, and its description in terms of quantum 
impurity problems. We consider a half filled Hubbard model on a fully frustrated 
lattice Q with a semicircular density of states in the limit of infinite lattice co- 
ordination. The term frustration refers to the degree of magnetic frustration in 
the parent Mott insulator. 
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Rcf Q reported that this model exhibits a Mott transition at a critical value 
of the ratio ¥ . The correct description of the destruction of the metallic state 
at zero temperature as the Mott transition is approached by increasing U, was 
proposed by Zhang, Rozenberg and Kotliar (ZRK) j^] on the basis of the iterated 
perturbation theory (IPT) S. 

The calculations of ZRK revealed that while the approach to the Mott tran- 
sition from the metallic side is driven by a collapse of the Fermi energy, as in 
Brinkman Rice theory, it also exhibits new unusual features. The metal 
disappears into an insulator with a preformed (finite) Mott Hubbard gap. We 
will refer to this as the semicontinuous scenario to be distinguished from the 
competing, bicontinous, scenario where the gap closes at the same point where 
the quasiparticle weight vanishes as discussed on page 0, 

The complete picture, of the Mott Hubbard transition in infinite dimensions 
emerged with the work of A. Georges and W. Krauth [jll| |l2| and Rozenberg 
et. al. juj, who described the destruction of the insulating state at zero 
temperature, the first order finite temperature metal insulator transition, and 
the crossovers that govern the behavior above the finite temperature second 
order critical point. They produced a wealth of physical results, which were in 
surprisingly good qualitative agreement with experimental data. [|l5| . The zero 
temperature scenario for the destruction of the metallic state was put on a firm 
footing by the development of the projective self consistent method ]l6| ]. This 
method overcame the difficulties associated with the presence of several energy 
scales, which had beset earlier treatments. 

In spite of these developments, several questions about the solution of the 
Hubbard model in large dimensions were raised [[l7| |l8| |19 p0[| a nd numer- 
ical studies were undertaken in an attempt to answer them. J21|. [ p2| [p3[ 
This renewed interest and in particular the insightful questions of Nozieres ]19| , 
motivates us to reexamine the problem from a new perspective, that of a Lan- 
dau -like free energy functional of a " metallic order parameter" , generalizing 
an approach used in our earlier studies of interacting random models with V. 
Dobrosavlevic |2^| . 

Our discussion highlights the peculiar character of the Landau theory of the 
metal to insulator transition. This singular dependence of the mean field free 
energy on the metallic order parameter (and not the specific approximations 
such as IPT, QMC or exact diagonalization of finite systems which were used 
in the early studies of this problem) is responsible for the unusual features of 
the solution of the Hubbard model in infinite dimensions. 

Our aim is partly pedagogical, we use the Landau functional to describe 
from a new perspective results that were obtained a few years ago. Besides 
clarifying the existing confusion in the literature of the subject, there is another 
purpose in writing a pedagogical note. There are not that many solvable models 
of the Mott transition in dimensions higher than one! We believe that there 
are still many lessons to be drawn from the solution in the limit of infinite 
dimensions, that can be of use in tackling more difficult problems, in the field of 
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strongly correlated electron systems. We believe that the Landau-like approach 
which we advocate in this paper can be valuable in other dynamical mean field 
studies. Finally, while we believe that the nature of the Mott transition in fully 
frustrated systems in the limit of infinite dimensions, has been understood at 
the qualitative level, there still remains a large amount of quantitative work 
to be done on this problem. Our insights should be a helpful guide to further 
investigations. 

This paper is organized as follows. After setting the notation in section 
||, summarizing the scenario describing the destruction of the metallic state in 
section ||. We state the questions raised by this suggestion and describe the 
alternative (bicontinous) scenario where the gap closes at the same point where 
the quasiparticle peak disappears in section |J. Two technical tools are essential 
to justify the validity of the semicontinuous scenario, the Landau free energy 
functional is described in section |[ and some results of the projective self con- 
sistent method are summarized in section ^. Using these tools, we describe the 
energetics of the metal insulator transition, inspired by the questions of Nozieres 
The Landau functional, provides us with a concrete bridge between the 
impurity model and the lattice model allowing us to use our knowledge of Kondo 
impurity physics to understand the Mott transition problem. 

In section [t| we use the Landau functional to describe the arguments of 
Fisher, Kotliar and Moeller |25| for the determination of the conditions for 
U c i, the point where the insulator disappears. Near U c \ the physical picture is 
that of an impurity in a weakly coupled regime, Nozieres exhaustion ideas are 
applicable in this case. 

In section H we recall the arguments of Moeller et. al. (T^| for the disappear- 
ance of the metal at the critical value U C 2- Here, the Mott insulator with a finite 
gap is indeed unstable towards the formation of a narrow metallic band at the 
Fermi level. The effective impurity description is in an intermediate coupling 
regime. From the perspective of our analysis based on a Landau functional, 
the semicontinous scenario, i.e. the fact that U c \ < U C 2, is an unavoidable con- 
sequence of the different behaviors of quantum impurity models in weak and 
strong coupling limits. 

In section |9| we argue that a more realistic consideration of the magnetic 
correlations in finite dimension, may change the character of the free energy 
functional and comment on the relevance of the dynamical mean field theory 
results to finite dimensional systems. 
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2 Lattice Model and Associated Impurity Hamil- 
tonian 



We consider the Hubbard model on the Bethe lattice in the paramagnetic phase 
with coordination d and hopping at half filling. 

<i,j>cr i 

The half bandwidth is given by D = It and we will use D = 1 as a unit of 
energy. The kinetic energy per site, K, can always be expressed in terms of the 
non local Green's function Gij. In the limit of large lattice coordination it can 
also be expressed in terms of the one particle Green's function: 

<K>=-±- ^GyW^^ff^GK (1) 

cr,<i,j>uj uj 

In the limit of large dimensions the total energy per site, E = < J ^ > reduces to: 
E = TJ2 [(«<> + li)G(iu) - l]e M0+ + i < K > (2) 

The interaction energy per site is given by: 

< U >= E- < K >= U < n^riii > (3) 

As is well known now jj| , all the local correlation functions of the model can be 
obtained from an Anderson impurity model with hybridization function 

am = y (4) 

provided that A(iuj) obeys the self-consistency condition: 

t 2 G(iu>)[A] = A(kj). (5) 

Here, G(mj)[A] is the Green's function of the SIAM (single impurity Anderson 
Model) 

e KC+ a C ka +J2 V * (4a U + U + Cka) + Uf + f ] f+f i = HsiAM- (6) 
Ka Ka 

viewed as a functional of the hybridization function A(ilu) which is the Hilbcrt 
transform of J2k ^k 2 5{ui — e&). Equations || and ^] express the total energy of 
the lattice model in terms of the local Green's function of the problem. We can 
therefore express the total energy in terms of the local spectral function p{u) = 
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-^ImG(iu> = u> + i5) using the spectral representation G(ioj n ) — J duj- f 
= J A^-e-niiu)) W ^ n ^( e ) tne semicircular lattice density of states: 



E = I f(w)(w + n)p{u>) + 2t 2 I I dutduiffa f^^ (7) 



UJ\ — L02 

We will work in the grand canonical ensemble with the chemical potential chosen 
to be equal to [i = U/2. f(u>) is the Fermi function. 



3 Evolution of the Spectral Function at Zero 
Temperature 

In this section, we describe the qualitative features of the evolution of the spec- 
tral function, as a function of interaction strength U/t which is obtained by 
solving the mean field equation at zero temperature. These features were 
discovered in an IPT || study by Zhang, Rozenberg and Kotliar. 

We start at large U with a paramagnetic insulating solution with a gap 
A(C/). When U is reduced below a critical value of U, denoted by U C 2 , (with 
A g = A(J7 c2 ) 0) the paramagnetic Mott insulator becomes unstable against 
the formation of a metallic resonance at zero frequency. 

The mathematical description of the ZRK scenario of the evolution of the 
spectral function when U C 2 is approached from below is the following: 

3.1 ImG(w, U) 7^ for all < and for all U < U C2 (finite spectral density 
everywhere in the metallic phase). 

3.2 \im u ^ u - lniG(uj,U) — for fixed u such that < \uj\ < -g- (Existence 
of a finite gap at the Mott transition point). 

We now discuss more delicate issues, in which the frequency approaches zero 
while at the same time, U approaches the critical value U C2 . More specifically, 
we define w = — — and take the limits w — > and u> — > such that x = to lw 
is fixed. This limit defines the scaling functions which were computed in ref 




3.3 linij^jj- Im G(wx, U) ^ 0, i.e. there is a finite density of states at the 
Fermi level all the way up to the transition. In particular the pinning con- 
dition which leaves the density of states at zero frequency unrenormalized 
is obeyed everywhere in the metallic phase. 

3.4 For a generic value of x lim,, „- ImEfe, U) is finite. 

Notice however that since Fermi-liquid theory is valid below the Fermi 
energy in the metallic phase, for a fixed value of U below U C 2 

lim JmS(cJ, U) = (8) 



G 



3.5 There exists aio~ 0(1) such that 



lim Im£(wo! U) = oo 



(9) 



This incipient divergence and its significance was recognized in ref |9[ : it repre- 
sents the precursors of the Hubbard bands in the metallic phase. Its presence is 
unavoidable, since spectral features resembling the Hubbard bands are already 
well formed on the metallic side of the transition jj| . The divergence of the self 
energy occurs outside the Fermi liquid regime and should not be interpreted 
in terms of quasiparticle scattering. It should be understood as the precursor 
of the pole found at zero frequency in the paramagnetic insulator phase. This 
pole indicates that the paramagnetic insulating ground sate, is not smoothly 
connected to the non interacting Fermi gas. 

It is important to stress, that in the metallic phase the density of states 
does not vanish for energies less than -S- . This is a simple consequence of the 
self consistency condition of the dynamical mean field theory. The statement 
that the Mott Hubbard gap is finite at the Mott transition point, should be 
understood in terms of the previously described, highly non uniform, limiting 
procedure. 

We stress that the results discussed above, were derived by non perturbative 
means. The mapping of the Hubbard model in large dimensions, onto the impu- 
rity model can be done using the cavity construction Q] which does not involve 
any expansion in U. Furthermore, to reach the conclusions discussed above, non 
perturbative treatments of the impurity model and the self consistency condition 
are required. 

In the next section we mention a perturbative expansion, the skeleton ex- 
pansion, which expresses the self-energy as a power series in U and in terms of 
the fully renormalized Green's function, 



Here, I a ^ m denotes the contribution of a specific Feynman skeleton graph, la- 
beled a and of order m in the interaction strength U, to the imaginary part of 
the self-energy evaluated at a frequency u>. 

The convergence properties of this series are not well understood |50| pT| . 
Since the Anderson impurity model with a hybridization function which is non 
vanishing at zero frequency, has a singlet ground state which is a smooth func- 
tion of U, it may converges for very small U. It is also known that the series 
diverges when U is sufficiently large and the lattice model supports a paramag- 
netic insulating phase. 

In the ZRK scenario, since the graphs of the the skeleton series for I a ,m(u) 
are evaluated in terms of U and G which has a very small spectral weight at 
low frequencies 




(10) 
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limy-^ ia.mM = ° (11) 

for all % > \w\ > 0, 

But Eqs. H H imply that the function which the skeleton expansion repre- 
sents in some form, behaves very differently in various frequency reanges. So 
even if the skeleton expansion converges pointwise in the open interval (0, U c i) 
the convergence in this interval cannot be uniform. Finally, we notice that ex- 
actly at the point U c i, the quasiparticle peak has zero weight. The system is in 
the paramagnetic insulating phase where the skeleton series is known diverge. 

The lack of uniformity in the frequency domain, is the mathematical mani- 
festation of the collapse of the Fermi energy, as we approach the transition. Be- 
low that scale a power series in the interaction has to be well-behaved because 
at low frequencies the system resembles a correlated metal, which is smoothly 
connected to the non interacting system by Fermi liquid theorems. At high 
frequencies, the system resembles a paramagnetic insulator, which has a doubly 
degenerate ground state at each site. For such a system skeleton perturbation 
theory is known to diverge, because a doublet cannot be smoothly connected to 
a singlet ground state. 



4 Critiques of the ZRK scenario 

Some of the findings in the ZRK paper described in the previous section were 
expected. For example, the gradual narrowing of the resonance as the Mott 
transition is approached, is the result of the Brinkman Rice mass enhancement 
[ p"o[ . Other aspects of the ZRK scenario, however, were new counterintuitive 
and surprising. The instability of a Mott insulator with a finite gap , towards 
metalization was unexpected (previously, such an instability was only expected 
to take place when the gap was infinitesimal). Also the incipient divergence 
of the self energy, at a relatively high energy scale y/ibD had not appeared in 
earlier slave boson studies. 

The alternative scenario for the Mott transition in frustrated systems, is a 
bicontinuous one (i.e. continuous from the metallic and the insulating side) In 
this scenario the gap closes from the insulating side at the same critical value of 
the interaction at which the resonance vanishes upon approaching the transition 
from the metallic side. This bicontinuous scenario was shown to occur within the 
slave boson formulation of Kotliar and Ruckenstein [^6| after including Gaussian 
fluctuations on top of the mean field theory p7| |28| and within a large N model 
of the metal to charge transfer insulator [|28| . The natural extension of this 
scenario to finite temperatures gives a smooth crossover between a metal and 
an insulator, excluding a first order phase transition between a metallic and an 
insulating phase but other extensions are possible [^9| . 

Several authors raised questions about the physical meaning and the internal 
consistency of the ZRK scenario and raised the possibility @ Q Q @ || 
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that the alternative, bicontinuos scenario might be realized in the frustrated 
large d Hubbard model with the semicircular density of states. 

Nozieres jl9| expressed the surprising nature of the evolution of the spectral 
function at zero temperatue in two different ways. First an energetic argument: 
if one is near the Mott transition, transfer of a small amount of spectral weight 
e from high to low energy across a Mott Hubbard gap A g costs eA g while the 
kinetic energy gain is only of order e 2 D. For the transition to take place, these 
two terms have to balance (i.e. A s has to be of order eD) thus ruling out the 
possibility of a preformed Mott Hubbard gap. 

The second argument makes use of Nozieres exhaustion principle ( which 
we interpret in the present context as saying that the Kondo quenching of a 
spin in a dilute bath of conduction electrons does not result in a substantial 
energy gain). Near the Mott transition, the ZRK spectral function has very few 
carriers, while there are many spins. The system is in the exhaustion regime, 
where very little Kondo energy can be gained, and is therefore unlikely to be 
energetically stable. Full Kondo compensation, which is necessary to form the 
resonance, and the Fermi liquid state is unlikely. 

Recently Kehrein argued [^0), that under the assumption that the sekelton 
expansion converges pointwise in the open interval (0, U e 2), the scenario pre- 
sented in the previous section is inconsistent. From equation [ll]hc concludes 

^ lim I a , m {uJ = x Vtii) = (12) 

He then interchanges the order of the infinte summation and the limiting pro- 
cedure. Using equation [l^ he obtains 

lmi u _ tU - ^(^/wxq, U) = contradicting eq. [| 

This exchange of limits, is not allowed, even when the skeleton expansion 
converges pointwise for all frequencies, and even if the total spectral weight of 
each graph can be bounded uniformly as is stated in footnote 14 of ref [E0| 



The exchange of limits is mathematically justified only when the series 10 
converges uniformly . The skeleton expansion cannot converge uniformly as 
one approaches the Mott transition since it diverges at U C 2, which in the ZRK 
picture, is in the paramagnetic insulating phase. 

Recently E Lange Q analyzed the soluble strong coupling limit of an An- 
derson impurity for arbitrary U. He showed explicitly that while the skeleton 
expansion may fail to converge, particularly, at high frequencies, an IPT-like 
treatment, gives the exact answer for this model. 

The rest of this paper addresses the questions of Nozieres. The goal is to 
explain why at some point it becomes energetically favorable to metallize a 
paramagnetic insulator while it still has a finite gap. We also clarify how the 
energetics of lattice models can be estimated using impurity models. For these 
purposes, we develop a framework based on a Landau like functional to connect 
lattice and impurity physics in the next section. 



5 Landau Functional 



To address the questions about energetics and to illuminate the analogies and 
differences between the dynamical mean field theory and the more standard 
Landau theory of phase transitions it is useful to introduce the free energy 
functional |2J : 

F LG [A] = -Tj2 + F imp [A] (13) 

u 

where Fi mp is the free energy of the impurity model defined by ^ which can be 
represented as a functional integral: 



Here, Li oc is the action of a local f level with the hybridization set to zero. 

One should regard equation |l3| as a Landau Ginzburg Functional of the 
"Metallic Order Parameter" A(ioj). At zero temperature, the function A(iui) is 
non zero in both phases, but it has very different low frequency behavior in the 
metallic and in the insulating phase. 

Differentiating the free energy with respect to A(iuj) we obtain the mean 
field equation [|. The derivative of the first term is — 2 A ^"- >T . Combining it 
with the results of differentiating the impurity free energy in [l4] 

rp 

mp =TS(T</(r +( icJ ) /(T ( iw ) >= 2TG(icj)[A] (15) 



5A(iuj) 



we obtaine eq. |[ 

We will show below, by means of an explicit calculation, that evalu- 
ated at the saddle point || gives the correct energy of the lattice problem. Flg[A] 
allows us to consider the free energy of arbitrary hybridization functions (i.e. 
away from self consistent solutions), and to interpolate between various station- 
ary points. 

It also has a nice physical interpretation by analogy to the mean field free 
energy of a spin system regarded as a function of the Weiss field h: 

0F La [h] =P^j~ log[ch[2(3h]] (16) 

The first terms in |l6| and 13 are the cost of forming the Weiss fields around 
a site. The second terms are the energies of a site (spin in the classical case, 
electron in the quantum case) in the the presence of the the Weiss field around 
it. 

Note that our order parameter A(iu) is always non zero but has a qualita- 
tively different form at low frequencies on either side of the transition. This is 
analogous to the density in the liquid gas transition. At zero temperature, there 
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is a qualitative difference between a metal and an insulator, which appears nat- 
urally in this formalism. The metallic phase is characterized by a hybridization 
function which is nonzero at low frequencies, as opposed to insulating phases 
for which the hybridization function vanishes at low energy |33 |. 

We stress the difference between the Landau functional [13 and expression 
0. Both [?| and [l3| give the correct energy if one evaluates them at the exact 
local spectral function or the exact hybridization function respectively. However, 
[l3] is a stationary functional, i.e. upon differentiation it produces the correct 
dynamical mean field equations. On the other hand, (7] is not stationary (its 
derivative with respect to p is nonzero at the physical spectral density) and 
therefore it does not give the correct mean field equations. It is not a Landau 
functional, it can only be used to obtain the ground state energy if one knows 
the correct spectral function but it cannot be used to determine the spectral 
function itself. 

We now prove, that the Landau free energy, evaluated at the saddle point 
value of the order parameter, indeed gives the total energy of the lattice model. 
In the process we make the correct connection between the various contribu- 
tions to the energy of the lattice model and to the energy of the corresponding 
impurity model. 

The lattice Hamiltonian contains two terms, K and U. When A(iuj) (or 
correspondingly the one particle Green's function) changes, both < K > and 

< U > change. Flc[A] has two terms, — TJ2uj A ^ an d F imp [A], but it would 
be incorrect to associate Fi mp [A] directly with the kinetic energy of the lattice 
model. At a stationary point and at T = 0, the correct connection is: 

F imp [A] = ^<K> + <U> (17) 

Notice that since we work at zero temperature, the free energy and the energy 
are identical. The impurity energy is composed of three parts, the local cor- 
relation < U > , the hybridization energy 2 J2ka ^ k < fl c <?k > which is given 
by 4T A(iu))G(iuj) (extra factor of 2 comes from spin) and the change in 
kinetic energy of the conduction electrons, which is defined as 

J2ka e k(< cl k C ak > A - < cl k C ak > A=0 ) 

and is therefore given by 



2T£ 



,k 



V k 2 e k 
(iuj - e k )' 



;G{iw) 



= -2T^A(wj)G(iw) + 2T^] )-[— -(A(iu)G(iuj)]iu = -T^A(iw)GH 

UJ UJ UJ 

Combining the three terms, we arrive at equation [l7]. The first term in the 
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functional in eq [D| is given by 

A (im) 2 -1 

1p~ ~ ~ 



(if) > (18) 



The positivity of this term is essential if we want to interpret it as the cost of 
forming the Weiss field. Clearly once the Weiss field is formed the impurity 
energy is negative. (Notice that we are in the half filled situation where A(io>)is 
purely imaginary). 

Notice that while local quantities have a straightforward interpretation in 
terms of the impurity model, nonlocal quantities such as the kinetic energy 
require more care. This is why we stated very clearly how the self-consistent 
local impurity approximation embeds the impurity model in a medium and gives 
well defined relations between the kinetic and potential energy of the lattice 



model, and the energy of the Anderson impurity model, eqs yj, [18| [13|, and [17 
are exact in the limit of large dimensions. 

The previous discussion highlights the self consistent character of the Local 
Impurity Self-Consistent Approximation, whereby the energy is composed of 
two terms, cost of the Weiss field and the energy gain of the impurity in the 
presence of the Weiss field. This is very different from the reasoning within the 
local impurity approximation, which would equate the energy of the lattice to 
the number of sites times the impurity energy. 

Nozieres Jl8| jl9| reduced the validity of the semicontinuos or the continous 
scenario to the following question: does one gain Kondo energy once per spin or 
once per electron in the resonance (described by the small fraction of electrons 
contained in the low energy part of the bath)? . 

To answer this question an unambiguous expression of the energy of the 
lattice models in terms of the energy of impurity models in the putative trial 
states is required. Expression |l3|, unambiguously describes the energy per lattice 
site, i.e. the energy per spin, in a self consistent local impurity framework 
substantiating the semicontinuous scenario. 

We can now describe in a precise manner the energetics near the critical 
points U c i and £/ c2 . Their existence highlights an essential difference between 
the classical and the quantum Landau Ginzburg functionals. While in the clas- 
sical case the Landau functional is an analytic function of the magnetization 
and of the Weiss field which has a straightforward power series expansion, this 
is not the case in the quantum case. 

Free energies obtained by "integrating out" massless degrees of freedom ( 
fermions with a Fermi surface) are non analytic. Here we will deal with the 
manifestation of this problem in infinite dimensions. We note however that 
while this problem is most extreme in infinite dimensions, some manifestation 
of this problem in weaker form may be relevant to the treatment of quantum 
phase transitions in finite dimensions. 

An important lesson to be learned from this analysis is that the quantum 
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Landau Ginzburg functional is non analytic in the field A. The source of non 
analyticity can be traced to the presence of a Fermi surface, and in infinite 
dimensions can be given a sharp formulation following Nozicrcs' analysis of the 
Kondo model. 



6 Projective Self Consistent Analysis 

In the next sections we address a number of questions concerning the stability 
of the paramagnetic insulating state A (ui). When we evaluate the free energy 
functional in the perturbed state A (v) + 5Al(ui) (i) do we gain or lose energy 
? (ii) and how much? (iii) What is the physical origin of the energy gain ? 

For this purpose we analyze the low energy behavior of the Landau func- 
tional, using the projective self consistent method, a technique designed to iso- 
late the low energy physics of impurity models in a bath which is self consistently 
determined (as in the solution of lattice models in the limit of large lattice co- 
ordination) 

The goal of this section is to set up the notation associated with this re- 
duction to low energies. The results quoted here will be used in the next few 
sections to isolate and discuss the singular dependence of the Landau free en- 
ergy on the metallic order parameter. The projective self consistent method 



was developed in references |16| ]25| and [|3J] where it was used to calculate the 
value of U C 2 and the scaling functions at the Mott transition point. 

We start with a Mott paramagnetic insulating solution, A (w), i.e. a solution 
to the mean field equations with a finite gap and we add to it a perturbation 
v l2 ■ 

SAl^ui) = ^2 ek ijL e localized in the low energy region (i.e. the variables e k 
are much smaller in absolute value than the Mott Hubbard gap). 

Notice that since Ao is the exact insulating solution, (i.e. a stationary point 
of the functional |l5| ) the expansion of the energy in the small addition 6Al(u>) 
starts with quadratic terms We will therefore only write the terms which are 
quadratic in SAl{uo), and ignore the linear variations which have to cancel in 
the final answer |p2[ . 

Furthermore we take SAl(uj) as entirely concentrated in the region of low 
energies, since changes at high energies will be shown to be unimportant (see 
discussion of this point on page [l6]). Using a Schrieffer Wolff ]35| transformation 
on the Hamiltonian, we can reduce ^ (to order V^ 2 ) to a Kondo like Hamiltonian 

H K (SA L ) =E ins + J s {U,t)Y^V^clac kl Vi!S + Y,t k ct (J c ka (19) 

Ei ns is the energy of the Mott insulator, the subscript k, runs over the low 
energy conduction electrons. We introduced the local conduction electron: 

- = I>fe (20) 
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and S, the local renormalized impurity spin S a — ^cr a aa /X aa i acting on the 
low energy spin degrees of freedom. | J.) , | |) a . These are defined formally as 
ground states | J,) , | |) , of the Anderson impurity model with hybridization 
function Ao X are Hubbard operators, acting on those spin degrees of 

freedom. 

J s (U,t) is a monotonically decreasing function of U/t which depends ex- 
plicitly on the insulating solution Aq. An explicit expression is given by: 

Since we intend to compute the energy change to order <5Al(zu;) 2 , we need 
in principle the effective Hamiltonian obtained by the Schrieffer Wolff ]35j[ trans- 
formation to order V fc ■ It has the form: 

<i= §A 3) S-S, L + ^EAcLc-eL, + 4 L(T C L<7 ) (21) 

+ § J 3 ^S- s L + ^J 4 (3) K T - \){n Li - |) (22) 
with c eL a = 2J2 Vk^kCka and s K = j(cl a a al3 cp + c^ a u a ^c^). The coefficients 

k 

have been evaluated numerically |m| but they are not necessary for our 
purposes since the expectation value of the fourth order Hamiltonian ^2] in the 
ground state is zero, hence contributions to the energy to order 8/S.i J (iuj) i only 
arise from a second order computation using eq. [l^. 

It is also instructive to record the explicit form of the low energy part of the 
Anderson model f operator after the canonical transformation is performed: 

F a = -^-DVzY,S a a a a,a>c La/ (23) 

a' 

Z is a very important quantity, the total low energy spectral weight of the 
trial state: 

^ = 4/^l(-')^' = E 4 C (24) 

J k 

It is proportional to the quasiparticle residue of Fermi liquid theory. The propor- 
tionality constant between the low energy spectral weight and the quasiparticle 
residue was evaluated in rcf |1| . 

We now isolate the dependence of the Landau functional on <5Ai(icj), from 
cq. |l3[ namely, 

E L (SA L ) = -I^,5A L M 2 + E K (6A L ) (25) 

UJ 

Ek is the usual expression for the energy of a Kondo impurity model 

E K {5A L ) =< H K > [V k L ]~ <H K > [V k L = 0] (26) 
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While one can make qualitative arguments about transfer of spectral weight 
by regarding El(SAl) as a function of Z one should always keep in mind that 
El is really a functional of the whole hybridization function A ( or, if we restrict 
ourselves to low energy variations, of 8Al). In fact, we shall show in the follow- 
ing sections that the free energy functional ( p^ ) takes very different values for 
hybridization functions having the same value of Z but different aspect ratios 
or shapes. The physical reason underlying this is the different behavior of the 
Kondo energy depending on whether the model |l^ is in the weak coupling or 
the strong coupling limit. This will have significant implications for the Mott 
transition in the limit of infinite dimensions. 



7 Energetics near U, 



cl 

In this section we describe the energetics of making a "very small " perturbation 
(very small is defined by the requirement that the corresponding effective Kondo 
problem |l9| is weakly coupled) to the Mott insulating solution. The goal is to 
determine the location of U c i, the point where the insulator ceases to exist 
because it becomes linearly unstable against such a small perturbation. p5| |. 
The cost of modifying the Weiss field (first term in ), is computed by 

v l2 

inserting explicitly the expression 8Ai,(iLo n ) = ^ fc ^ k _ ek , doing the Matsubara 
sums, and taking the zero temperature limit, and is given by: 

o v l 2 v l 2 

e k ,>o,e k <0 

Next we evaluate eq. ^6| in second order perturbation theory to compute the 
energy gain 



E K (5A L ) = - Y ( JsV k V v) 2< fo\S-4*vck'*c+ a ac k J\fo > 

|$o > is a Fermi sea, i.e. a ground state of J2k £ kC + kaCka- Evaluating |28| we get 

E K (6A L ) = J s 2 |E £l . i>oe k <o (e.,-el) Therefore, the energy gained by a weakly 
coupled Kondo impurity embedded in the perturbed Weiss field is 

E K (SA L ) = -T[J]<5A L M 2 ]|j J s 2 (29) 

Balancing cost and gain then gives the condition for U c i, namely ^ = | J 2 (U, t). 
Recalling that t = ^ we have rederived the equation that determines U c x first 
obtained in ref |E5|: 
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5* = J ^r 6 (30) 

Since the closure of the gap, necessarily implies instability towards metal- 
ization, the U that solves eq |3(] is an upper bound to the value of U at which 
the gap closes |3g| . 

An advantage of the Landau formulation is that the trial states 5Al(zo;) 
do not have to be normalized to unity i.e. the condition 1 = ^ w A la e 1 " 
does not have to be satisfied. However, if for esthetic reasons one wants to 
exhibit trial states which do satisfy a proper normalization, this can be easily 
done without altering significantly the energy of the state by adding a small 
perturbation 5Ah centered at high energies near the gap edge. Since Ao was 
stationary, the cost associated with this addition is of the form 

a£ ,KS.,^W„W (3l) 

with a\ and 02 numbers of order unity. Since Vh 2 is of order Z, but the 
denominator D is much larger than the corresponding denominator in equations 
p7| and |2§| it is clear that this term does not contribute substantially to the energy 
balance and the condition for the insulator to become unstable is still controlled 
by the equation involving low energy denominators (^, p8| ) . 

The first argument of Logan and Nozieres, which stated that we can only 
gain energy from the metalization process when the gap has closed is in com- 
plete agreement with our detailed calculations, but is only valid for trial states, 
relevant to the determination of U c i. Under the assumption that the trial state 
is in the weak coupling regime, we can regard the changes in the spectral func- 
tion as small continuous deformations of the insulating solution. This kind of 
perturbations destabilize the insulating solution only when the gap closes. 

We can now address the second point raised by Logan and Nozieres. What is 
going on from the point of view of Nozieres exhaustion principle: a few electrons 
cannot gain a great deal of Kondo energy when they have to screen a lattice 
of spins. We see that this argument is qualitatively correct when applied to 
the dynamical mean field theory provided that the conduction band is weakly 
coupled to the spin. Indeed, if we regard the variables in eqs [27] and ^9] to 
be of order D, (and not of order ZD), the changes in energy are indeed of order 
Z 2 D 

More generally, we note that the partial screening of magnetic moments, is 
a very common situation in the full solution of the non linear dynamical mean 
field equations ^. It occurs in the correlated regime at finite temperatures when 
the condition A(iO + ) = — it is strongly violated. 
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8 Energetics near U C 2 



A most interesting aspect of the large d solution is the instability of the Mott 
insulating state with a finite gap to the formation of a narrow metallic band 
at the Fermi energy. This takes place at a critical value U C 2 > U c \, as we now 
show. 

We are interested in the energetics of this problem which we shall consider 
from the point of view of the free energy |?j and [l3|. We first consider the 
energetics of the problem by inserting a spectral function p(u>) = pa(u>) + 8p(uj) 
in the functional ^. po describes the Mott insulating state, and Sp is comprised 
of two pieces, one removing spectral weight Z from the gap edge at — A g and 
one that adds it at zero energy. 

Since the functional [7] is quadratic, the evaluation of the energy difference 
between this state (our candidate for metallic state) and the Mott insulating 
state po is straightforward. It consist of a linear part in Sp, SEi and a quadratic 
part in Sp, SE2 . These are given respectively by: 

SE 1 = f du 1 8p(u 1 )[u 1 f(Lo 1 ) + 2t 2 [ dw 2 p (u 2 ) f ^ ~ (32) 
SE 2 =2t 2 [ [ du*d^fW Sp ^ )Sp{U2) (33) 

J J LUl — LU 2 

It is now straightforward to evaluate SEi , which as expected turns out to be 
positive representing a cost of energy and given by 

r°° —A 

8E 1 = Z[A g + 2t 2 de Po (e) - 9 ] (34) 

JO [ e + ^£fJ e 

We can now see the breakdown of the exhaustion arguments when applied 
to the impurity model relevant to U C 2 . If Sp were indeed a "small" (in the 
sense of the previous section) perturbation, it would then be true that the 
contribution from SE2 is negative but small, of the order Z 2 and that to obtain 
a net gain from 5E\ +SE2 one would require the gap to be of order Z rather than 
finite. However the perturbation that we are considering (which is a perfectly 
legitimate trial state, for the large d Hubbard model) is not small. The low 
energy part (Sp = SpL + Spu) has height one and is described by the scaling 
function 8pl{^>) = j> c f'(~£j))- The contribution from 8E2 is of order Z and is 
given by 

5E 2 =Zt[ [ d^fi ^ 1 ^ (35) 

J J LJl — Ld2 

Now it is completely clear that the gain in kinetic energy, 8E2 is of the same 
order as the positive contribution 8E\ so an instability to the metallic state 
when the insulator has a finite gap can take place. To see that it does take place 
we need to resort to the variational functional O. 
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To demonstrate this point we take an insulating solution with a finite gap 
with U slightly bigger than U c i- According to the analysis of the previous 
section this is stable against a perturbation SAl which is in the weak coupling 
regime. On the other hand this state can be unstable against a perturbation 
SAl described by a Kondo model in an intermediate coupling regime. In this 
section we will concern ourselves with a perturbation 5Al(u>) with height of 
order unity . 

We will take the perturbation with the scaling form 

SA L (u)=DM-^) (36) 

with </>i(0) of order unity. It is instructive, following reference ^M , to write 
the low energy Kondo Hamiltonian corresponding to such a perturbation in a 
way that makes completely clear that the condition of unit height of SAl(Q) 
corresponds to an effective impurity model in the intermediate coupling regime 
(Kondo energy comparable to electron kinetic energy). 

H K = ZJ s (U,t)c L +ac L S + Y, £ kct a c ka (37) 

o -ZD<e k <ZD 

It describes a band of narrow electrons with bandwidth ZD interacting with an 
orbital the local electron of equation which is properly normalized ({c^cl} = 
1), via a Kondo exchange ZJ s (U,t). The Kondo problem in ^ has a coupling 
J e ff — ZJ s (jy) and an effective bandwidth D e ff = ZD. As a result "grjj 
is independent of Z, so the Kondo impurity is in the intermediate coupling 
regime. Therefore (SlS) < where we defined the local conduction electron 
spin operator by 

Sl = \c) L ac L (38) 

Notice that, while the value of the expectation value of the scalar product of 
the impurity spin and the local spin operator depends on the full form of the 
scaling function and has only been calculated numerically the fact that it 
is non zero depends only on the intermediate coupling nature of the associated 
impurity model in eq |37]. 



We now insert the variation Atrial = Ao + SAl, with SAl given by Eq. 36 
into the Landau function ( |l3| ) and compute explicitly the quantity 

m^lim ^ £G j At " a ' ] (39) 

Z^Q dZ ^ ' 

If m is positive, there is a finite net cost to forming the resonance, if it is 
negative then it is energetically favorable to metallize. From the chain rule, we 
find 

m = - 2T Y f^ Wn) ~ G L {iu n )}±5A L {iu n ) (40) 
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We now insert the spectral representation of the low energy part of the 

(iw n -e) ' 



impurity Green's function Gl(«^)[<5Al] = / de ^^^ , with pz(e) = ^4>(yd), 
into cq. ^ to obtain 



' v 1 r ,D, 



m = AD I dxj dyj^—A{jYct> l {x)-ct>[x)] ( t> x {y)y (41) 



o 



We next assume a trial state such that (j>\ (y) = at small y and (j>\ (y) = c 
between Wl and W2. Then eq. |4l] simplifies to: 



m 



Wl ^ 



W2 



4Dc dyy dx- r[40i(x) - <f>(x)] (42) 



(x - y) 



For Wl and W2 large, we thus arrive at 



m 



AcD(Wl — W2) [ dx[A<t>x(x) - 4>{x)\ (43) 

J — oo 



4> is of course a complicated functional of <p\. However, to estimate a lower 
bound to the instability point we only need the integral of <fi which we calculated 
in appendix B. Using those results, and the fact that both <\>\ and 4> are even 
we obtain the condition for U C 2 first derived in ref |l6| ]: 

1 -i5^! [ |_ < M I>1< (44, 

Here, Sl = \c + Lo& a a,a'CLa' is the local conduction electron spin. 

The proof of the instability of an insulator with a finite gap results from 
comparing conditions |3^ and ^ which imply the strict inequality: 

U cl < U c2 (45) 

This follows once we recognize that while (S.Sl) depends on and can 
only be computed numerically after solving a non linear self consistent problem, 
the fact that it is non zero and negative is immediate, from the fact that it 
describes an impurity in the intermediate coupling regime. Since J S (U) is a 
monotonically decreasing function of increasing U, the solution of eq pffl has a 
lower value than the solution of equation 

Notice that the derivation is free from ad-hoc assumptions. The projective 
self consistent analysis is applied in the region where this method is valid, i.e. 
for a bath with clear separation of energy scales. The formulation in terms of 
the Landau like functional ^ is useful, because it allows the introduction of 
trial states, which can be analyzed more easily than the solution of the full non 
linear dynamical mean field equations ^. 

We end the section by recapitulating the logic leading to the inequality [45[ 
We start with the Mott insulator and slowly decrease the value of U. In section 
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we asked the question when is the insulating solution unstable to a small 
perturbation. The answer is an implicit equation for U, eq. Its solution, 
which we denote U c i is certainly bigger or equal than the value of U at which 
the gap closes, since when the gap is closing the insulator is definitely unstable 
against small perturbations. Then we consider a trial state corresponding to an 
impurity model in an intermediate coupling regime. The condition for instability 
against this specific trial state, is given by eq fl4| While the specific value of 
< S.Sl > depends on the trial state, the fact that the Mott insulator is unstable 
at a value of U strictly larger than Ucl and therefore posessing a finite gap, relies 
only on the fact that < S.Sl >< 0. This is always true for impurity models 
which are in the intermediate coupling regime. The trial state we used gives 
us a lower bound to the true value of U C 2, since improving the scaling function 
might result in a better trial state which can destabilize an insulator at an even 
larger value of U. (This optimization was done numerically in rcf flti)] ). 

We therefore have proved that an upper bound for the U which is needed 
for gap closure is strictly smaller than a lower bound for the U at which the 
insulator becomes unstable against metalization. 

Within the context of the Landau analysis we do not have to worry about 
the normalization of 8A.l(lo). To exhibit the instability against a normalized 
trial state, one can proceed as we did on page [l6|. 

Notice the two significant differences with respect to the previous section: 

a) The cost of forming the required Weiss field (while still given by the first 
term in eq. ^ and therefore having an identical form to that considered in the 
previous section ( see eq. p7j)) is now of order Z (and not Z 2 as in the previous 
section ) because the height of the spectral function at the origin is of order 
unity. 

b) The gain in energy from the free energy of the impurity model is now also 
of order Z This is not surprising because each individual term in the Hamiltonian 
of the Kondo impurity is of order Z. 

The difference between the determination of U c \ and U C 2 can then be traced 
to the different behavior of the impurity model in the weak coupling and the 
intermediate coupling regime. To illustrate this idea in a much simpler setting, 
we analyze a very simple free energy function in Appendix A. 

9 Implications and Outlook for Finite Dimen- 
sions 

Our analysis so far was confined to zero temperature. However the proof of 
existence of two coexistent solutions (one metallic and one insulating) in a finite 
interval of interaction strength at zero temperature, has immediate consequences 
for the finite temperature phase diagram. Away from the phase boundaries, the 
Landau functional, is a smooth functional of the hybridization function. Since 
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we have two coexistent solutions, at a finite distance of each other J37[ , at zero 
temperature, by continuity the two solutions persist at finite temperature. 

At finite temperatures there is no qualitative difference between the metal- 
like and the insulating-like solution. Just the low frequency density of states is 
quantitatively different (large in the metal, small in the insulator). The coexis- 
tence of two solutions at zero temperature results in a first order phase transition 
at finite temperatures since entropy favors the paramagnetic insulating solution. 

In the temperature interaction strength phase diagram of the frustrated 
Hubbard model refs |Q (l^] Jl3| Q the points U c \ and U C 2 are the end points 
of two lines U c \ (T) and U C 2 (T) which cross at a finite temperature second order 
point at {Umit,T M it)- 

Rozenberg et. al. |ll| observed that the U c i(T) and U C 2(T) can be contin- 
ued above Tmit, where they become crossover lines. The physical interpretation 
[ jl5| , as crossover lines, was remarkable similar, to that of the zero temperature 
points where they originated. At U C \(T) the gap becomes comparable to the 
temperature and activation results. At U C 2(T) the Fermi liquid coherence dis- 
appears 

The limit of full frustration is rare. However, it turns out that for reasons 
having to do with orbital degeneracy, crystal structure, and longer range multi- 
spin interactions, many three dimensional transition metal oxides undergoing a 
metal to insulator transition, are very nearly frustrated, and finite temperatures 
easily stabilize the paramagnetic insulator phase. 

The continuation of the zero temperature paramagnetic metal to param- 
agnetic insulator transition to finite temperatures, then becomes immediately 
relevant for answering a long standing question of whether one can account for 
the phase diagram of NiSe2- x S x or of V2O3 in purely electronic terms without 
including explicitly the coupling to the lattice. The answer suggested in refs [[ll| 
15 was positive, provided one includes a sufficient degree of magnetic 
frustration. 

The original discussion relied on approximate methods such as the IPT, 
QMC, or exact diagonalization, which were open to criticism. We have stressed 
in this article that at zero temperature there are two coexisting solutions in a 
range of interaction strengths U, which allows us to conclude (without recourse 
to approximations), that the qualitative features found in references fl5|| in 
IPT are not artifacts of the approximation, but genuine features of the solution 
of the dynamical mean field equations The character of the phase diagram, is 
a consequence of the analytic structure of the free energy functional. One can 
view, the different approximation methods, as providing different approximate 
values to the coefficients of the Landau Ginzburg functional of the order param- 
eter. As long as the functional has the correct analytic structure, the qualitative 
features of the solution are preserved (of course, approximate methods can only 
give approximate values of Tmit, U c i(T) U C 2(T) etc ). Approximations such as 
the fully self consistent skeleton method fail because they do not capture the 
correct analytic properties of the free energy functional, i. e. they miss the 
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paramagnetic insulating phase. 

The dynamical mean field solution of the Hubbard model has allowed us to 
study and understand the finite temperature consequences of the existence of a 
Brinkman Rice quantum critical point JlO[ ] U C 2, where spatial coherence is lost. 
This quantum critical points is always inaccessible because no system is fully 
frustrated and some form of magnetic order always sets in before reaching the 
critical point at zero temperature. The DMFT has been of extraordinary value 
in revealing the hidden origin of the anomalous temperature dependence [[l5| 
of many properties. The observable anomalous finite temperature properties, of 
the paramagnetic phases close to the Mott transition point, are connected, via 
DMFT, to the presence of a zero temperature quantum critical point. 

Having established the physical relevance of the paramagnetic metal to para- 
magnetic insulator transition in magnetically frustrated systems, we turn to the 
broader and more interesting question of how reliable the DMFT description of 
this transition is. 

Conventional mean field theories of classical phase transitions are qualita- 
tively valid from infinite dimensions down to a finite upper critical dimension. 
The situation is very different in dynamical mean field theory. As currently for- 
mulated, in the paramagnetic insulating phase, this theory omits a physically 
relevant scale, the magnetic exchange energy. The dynamical mean field picture 
of metalization, as a result, cannot capture the competition between intersite 
magnetic correlations and the Kondo effect. 

In a simplistic picture of the DMFT, which views the central site as a strongly 
correlated electron f, and the neighbor as a non interacting bath electron c, the 
f site will have a well formed moment S ( which captures the localized charac- 
ter of the electron) while the nearby sites are represented by a conduction band 
which is narrowed to represent the correlated but itinerant character of the orig- 
inal electrons. Since the bath is uncorrelated, it can only describe a moment, 
by becoming infinitely narrow. This does not happen in reality, as a result of 
magnetic correlations, which in three dimensions, ultimately trigger some form 
of magnetic order. The qualitative success of the dynamical mean field ap- 
proach depends crucially on the weakness of the intersite magnetic correlations 
compared to the renormalized Fermi energy tp *. 

Since in the dynamical mean field theory, the effective Fermi energy van- 
ishes at zero temperature and at the critical point, the situation at first sight 
looks hopeless. However the experimental confirmations of the existence of an 
underlying quantum critical point at U C 2 arise from the finite temperature be- 
havior. As long as the degree of frustration and the temperature is sufficiently 
high so that the relevant magnetic energy scale is smaller than the temperature, 
(or when the distance from the critical point is sufficiently large so that the 
renormalized Fermi energy is larger than the magnetic exchange), we are in a 
situation where many of the qualitative predictions of the dynamical mean field 
theory are applicable, some of these predictions have actually been observed 
experimentally |f5[ |38| • Underlying these successes are some unique features of 
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the dynamical mean field theory which are not shared by ordinary mean field 
theories. A single dynamical mean field solution can contain information about 
a large number of energy scales. As the corresponding dynamical mean field 
theory of spin glasses, it works well in frustrated situations. This should be 
contrasted with static mean field theories, which usually do not work well when 
there are several competing ground states. 

We now turn to the outstanding open issue of how the findings of the dynam- 
ical mean field theory are modified when one is at sufficiently low temperatures, 
or sufficiently close to the transition so that the magnetic correlations have 
to be taken into account. The considerations presented in this note suggests 
significant departures from the results of ynamical mean field theory. 

In infinite dimensions, near U c i we have complete decorrclation of S and 
Sl = cjjoch. Near XJ c i however we have a finite value of < S.Sl > because 
a band which is infinitely narrow, can be polarized without any cost, and we 
succeed in getting some finite energy from the spin correlations between neigh- 
boring sites all the way to the transition. 

This can not be accomplished once we have a finite exchange among nearby 
sites (an effect which is formally l/d and therefore invisible in the large d para- 
magnetic phase but becomes of order one once we allow for magnetic long range 
order) . Now the Kondo interaction between the central spin S and Sl will 
compete with the exchange interactions with the nearby spins. As a result the 
free energy becomes a smoother function of the hybridization. If there is some 
kind of Kondo effect going on at the transition it will now occur in a regime of 
weaker coupling strength. The same conclusion can also be reached from the 
point of view of the itinerant electrons. If we now accept that the self energy can 
contain k dependence the pinning condition p(Q) — is now eliminated. The 
possible Weiss fields can now be more like those considered in the discussion in 
section 0. 

For these reasons we believe that the destruction of the insulator near U c i 
which we described in section ^ may play a more fundamental role in finite 
dimensions or in partially frustrated situations. 

Stated in a different language the large d picture is one where ep = Zt, 
the renormalized Fermi energy, is larger than J^, the magnetic superexchange 
between a pair of spins. Clearly this is justified in the formal infinite d limit 
because = tfj/U oc ^ is vanishingly small no matter how small ep is. That 
situation changes immediately once Jij is finite. Clearly new physics is expected 
when is larger than ep. 

Within DMFT the influence of magnetic correlations on the Mott transition 
has been recently studied recently by approaching the metal insulator tran- 
sition from an ordered state. In this case, one can quantitatively study how the 
magnetic correlations dramatically modify the behavior near the metal insulator 
transition point. In the absence of magnetic long range order, the task is more 
difficult. We have to do it without the guidance of the infinite d limit , because 
all those effects disappear in that case. These effects appear next order in l/d, 



23 



but a complete treatment of those corrections, is still missing 




We believe that the basic idea of this paper, i.e. that different types of Mott 
transition should be viewed as different types of bifurcations, of a system of 
functional equations for a metallic order parameter closely related to the one 
particle Green's function, will also be useful beyond the limit of large lattice 
coordination. 



10 Appendix A 

To illustrate the main features of the Landau Functional we discuss in this 
Appendix a toy free energy which has a structure similar to eq. [l3| but which 
is much easier to analyze. This free energy is a function of two variables : j 
(which should be thought of as one of the V k 2 of the previous sections) and t 
which should be thought of as a single e&. ft is given by: 

F[t, j] = ^ + {t- ^f.fJs 2 +t*) - (46) 

The analogy with the discussion in previous sections should be apparent, the 
first term is the cost of making the Weiss field, the second term is the energy 
that the impurity gains once the Weiss field is formed . We have mimicked 
the non analytic dependence of the free energy on the coupling strength in 
weak coupling by the exponential term which does not have a Taylor series in j , 
around j=0 The analysis leading to the determination of U c \ consists of simply 
evaluating whether there is an energy gain by increasing j away from zero, in 
weak coupling, that is for a fixed value of t and for an infinitesimal value of j. 
This stability analysis is equivalent to checking whether the second derivative of 
the free energy with respect to j, at fixed t is positive or negative at the origin 
From equation ^ we conclude that the U c i like instability occurs at the value 
of U ( J s is a function of U ) such that 

2 I 2 

--±= = (47) 

We now perform the determination of the critical value which parallels Uc2, i.e. 
we perform a minimization of the free energy without assuming whether the 
minimum occurs at weak or strong coupling. In other words, we focus on the 
functional which tells us if energy is gained or lost by increasing w , the weight 
at the origin, for fixed value of the inverse coupling strength a which we define 
as a = |. U C 2 is determined by then optimizing with respect to the coupling 
strength a . 



n ■ , l r i w 10J S W 

j[j = w,t — w.a\ = g[w,a\ = haw y ar wr + J s w z (48) 
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Now we notice that since in this toy model the free energy is linear in w, the 
instability occurs when this term is negative (w by definition is always positive) . 
The calculation that parallels U C 2 requires us to minimize with respect to a and 
search for instability. Terms non linear in w could be added to the toy free 
energy so as to obtain well defined solutions to the full non linear problem, but 
will not be considered here. 

It is clear that at U c %, the system prefers to be metallic and favors a strong 
coupling state. 



11 Appendix B 

To derive eq. ( fii]) we need the asymptotic behavior of the real part of the 
Green's function in the Anderson impurity model in the regime ZD <C cj <C A 9 , 
(i.e. for energies much less than the Mott Hubbard gap but greater than the 
resonance width ). From the asymptotic behavior 

5A L (iu> n ) « (49) 

it follows that f°° (j>i(x)dx = | The low energy part of the f electron operator 
is given in equation |23|. Using the definition of the local electron spin ^ and 
the identity 

i V < {(ScW, (SSc L ) a ) >= 2(|- < S.S L >) (50) 

a 

we obtain 

J p L (e)de =\J2< ( F ° F *°} >w f (^(g- < S-Sl >) (51) 
(j 

It follows that 4>(x)dx = ^(DJ S ) 2 (|— < S.Sl >) which is then used in eq. 
H| in the text. 

NOTE ADDED: After this paper was completed, two numerical studies of 
the problem discussed in this note have appeared. In an algorithmic tour de 
force, Bulla obtained a value of U C 2 in excellent agreement with earlier 
results p6| . Using a different numerical method Gebhard and Noack pi], ob- 



1. 



taincd a value of U c \ which is in excellent agreement with earlier work 1 15 
Their value for U C 2 however, is much smaller than all previous estimates. 
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